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We recently introduced the concept of "fc-essence" as a 
dynamical solution for explaining naturally why the universe 
has entered an epoch of accelerated expansion at a late stage 
of its evolution. The solution avoids fine-tuning of parame- 
ters and anthropic arguments. Instead, fc-essence is based on 
the idea of a dynamical attractor solution which causes it to 
act as a cosmological constant only at the onset of matter- 
domination. Consequently, fc-essence overtakes the matter 
density and induces cosmic acceleration at about the present 
epoch. In this paper, we present the basic theory of fc-essence 
and dynamical attractors based on evolving scalar fields with 
non-linear kinetic energy terms in the action. We present 
guidelines for constructing concrete examples and show that 
there are two classes of solutions, one in which cosmic accel- 
eration continues forever and one in which the acceleration 
has finite duration. 



I. INTRODUCTION 

A concordance of cosmological observations^ of large- 
scale structure, the cosmic microwave background 
anisotropy and Type IA supernovae at deep red shift 
suggest that the matter density of the universe comprises 
about one-third of the critical value expected for a fiat 
universe. The missing two-thirds is due to an exotic dark 
energy component with negative pressure that causes the 
Hubble expansion to accelerate today. One candidate for 
such a component is a cosmological constant (A) or vac- 
uum density. Another possibility is a dynamical compo- 
nent whose energy density and spatial distribution evolve 
with time, as is the case for quintessence^ or, as explored 
herein, for fc-essence. til 

A key challenge for theoretical physics is to address 
the cosmic coincidence problem: why does the dark en- 
ergy component have a tiny energy density (C(meV )) 
compared to the naive expectation based on quantum 
field theory and why does cosmic acceleration begin at 
such a late stage in the evolution of the universe. Most 
dark energy candidates (such as the cosmological con- 
stant) require extraordinary fine-tuning of the initial en- 
ergy density to a value 100 orders of magnitude or more 
smaller than the initial matter-energy density. Propo- 
nents of anthropic modelsQ often pose the problem as: 



why should the acceleration begin shortly after structure 
forms in the universe and sentient beings evolve? If the 
dark energy component consists of vacuum density (A) 
or quintessences in the forms that have been discussed 
in the literature to date, the answer is either pure coin- 
cidence or the anthropic principle. 

The purpose of introducing fc-essence is to provide a 
dynamical explanation which does not require the fine- 
tuning of initial conditions or mass parameters and which 
is decidedly non-anthropic. In this scenario, cosmic ac- 
celeration and human evolution are related because both 
phenomena are linked to the onset of matter-domination. 
The fc-essence component has the property that it only 
behaves as a negative pressure component after matter- 
radiation equality, so that it can only overtake the matter 
density and induce cosmic acceleration after the matter 
has dominated the universe for some period, at about the 
present epoch. And, of course, human evolution is linked 
to matter-domination because the formation of planets, 
stars, galaxies and large-scale structure only occurs dur- 
ing this period. A further property of fc-essence is that, 
because of the dynamical attractor behavior, cosmic evo- 
lution is insensitive to initial conditions. 

The existence of attractor solutions is reminiscent of 
quintessence, models based on evolving scalar fields with 
exponential "tracker"BB potentials. In these models, 
an attractor solution causes the energy density in the 
scalar field to track the equation-of-state of the domi- 
nant energy component, be it radiation or matter. An 
advantage is that the cosmic evolution is insensitive to 
the initial energy density of the quintessence field, and, 
for many models, the scenario can begin with the most 
natural possibility, equipartition initial conditions. (For 
the case of vacuum energy or cosmological constant, the 
vacuum energy must be set 120 orders of magnitude less 
than the initial matter-radiation density.) However, so 
long as the field tracks any equation-of-state, it cannot 
overtake the matter-density and induce cosmic acceler- 
ation. Indeed, for a purely exponential potential, the 
field never overtakes the matter density and dominates 
the universe. Hence, this is an unacceptable candidate 
for the dark energy component. In tracker models, the 
problem is addressed because the curvature of the po- 
tential ultimately dips to a critically small value once 
the field passes a particular value, Q such that the field 
Q becomes frozen and begins to act like a cosmological 
constant. The value of the potential energy density at 



1 



Q = Q determines when quintessence overtakes the mat- 
ter density and cosmic acceleration begins. The overall 
scale of the potential must be finely adjusted in order 
for the component to overtake the matter-density at the 
present epoch. So, while tracker models allow equiparti- 
tion initial conditions, they require the same fine-tuning 
as models with cosmological constant. 

The distinctive feature of the fc-essence models we con- 
sider is that fc-essence only tracks the equation-of-state 
of the background during the radiation-dominated epoch. 
A tracking solution during the matter-dominated epoch 
is physically forbidden. Instead, at the onset of matter- 
domination, the fc-essence field energy density e drops 
several orders of magnitude as the field approaches a 
new attractor solution in which it acts as a cosmologi- 
cal constant with pressure p approximately equal to —s. 
That is, the equation-of-state, w = p/e, is nearly -1. 
The fc-essence energy density catches up and overtakes 
the matter-density, typically several billions of years af- 
ter matter-domination, driving the universe into a period 
of cosmic acceleration. As it overtakes the energy density 
of the universe, it begins to approach yet another attrac- 
tor solution which, depending on details, may correspond 
to an accelerating universe with w < -1/3 or a deceler- 
ating or even dust-like solution with —1/3 < w < 0. 
In this scenario, we observe cosmic acceleration today 
because the time for human evolution and the time for 
fc-essence to overtake the matter density are both sever- 
als of billions of years due to independent but predictive 
dynamical reasons. 

The fc-essence models which we have found rely on 
dynamical attractor properties of scalar fields with non- 
linear kinetic energy terms in the action, models which 
are unfamiliar to most particle physicists and cosmol- 
ogists. Some of the concepts were first introduced to 
develop an alternative inflationary model known as fc- 
inflation.tj In this paper, we present a thorough, pedagog- 
ical study of dynamical attractor behavior and the appli- 
cation to present-day cosmic acceleration. The paper is 
organized as follows: In Section II we derive the basic 
equations describing the dynamics of a universe filled by 
matter, radiation and fc-essence. In Section III, we clas- 
sify the possible attractor solutions for fc-essence. In some 
cases, the attractor solution causes fc-essence to mimic 
the equation-of-state of the background energy density; 
we refer to this as a tracker solution. In other cases, 
fc-essence mimics a cosmological constant, quintessence 
or dust without depending on the presence of any ad- 
ditional cosmic energy density. In Section IV, we show 
how these principles can be used to control how fc-essence 
travels through a series of attractor solutions as the uni- 
verse evolves beginning from general initial conditions. 
In particular, we show how fc-essence can transform au- 
tomatically into an effective cosmological constant at the 
onset of matter-domination, as is desired to explain nat- 
urally the present-day cosmic acceleration. In Section V, 



we show how to utilize these concepts to design model 
Lagrangians. We explore two illustrative examples. In 
one case, the future evolution of fc-essence causes the uni- 
verse to accelerate forever. In the other case, fc-essence 
ultimately approaches an equation-of-state correspond- 
ing to pressureless dust, and the universe returns to a 
decelerating phase. 



II. BASICS OF i^-ESSENCE 

The attractor behavior required for avoiding the cos- 
mic coincidence problem can be obtained in models with 
non-standard (non-linear) kinetic energy terms. In string 
and supergravity theories, non-standard kinetic terms 
appear generically in the effective action describing the 
massless scalar degrees of freedom. Normally, the non- 
linear terms are ignored because they are presumed to 
be small and irrelevant. This is a reasonable expecta- 
tion since the Hubble expansion damps the kinetic en- 
ergy density over time. However, one case in which the 
non-linear terms cannot be ignored is if there is an attrac- 
tor solution which forces the non-linear terms to remain 
non-negligible. This is precisely what is being considered 
here. Hence, we wish to emphasize that fc-essence mod- 
els are constructed from building blocks that are common 
to most quantum field theories and, then, utilize dynam- 
ical attractor behavior (that often arises in models with 
non-linear kinetic energy) to produce novel cosmological 
models. 

Restricting our attention to a single field, the action 
generically may be expressed (perhaps after conformal 
transformation and field redefinition) as 
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where we use units such that =£k = l ; and 



x = -(y<py 



(i) 



(2) 



The Lagrangian p depends on the specific particle the- 
ory model. In this paper, we consider only factorizable 
Lagrangians of the form 



p = K(ip)p(X), 



(3) 



where we assume that K ((f) > 0. 

Lagrangians of this type are general enough to ac- 
commodate slow-roll, power-law and pole-like inflation, 
and they also appear rather naturally in the effective 
action of string theory. For small X, one can have 
p(X) = const. + X + 0(X 2 ). Ignoring quadratic and 
higher order terms, the theory corresponds (after field 
redefinition) to an ordinary scalar field with some po- 
tential. Normally, higher order kinetic energy terms are 
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ignored under the assumption that they are small, but 
the attractor solutions considered here insure that the 
non-linear terms remain non-negligible throughout cos- 
mic history. The scalar field for which these higher order 
kinetic terms play an essential role we call, for brevity, 
k— essence. 

To describe the behavior of the scalar field it is con- 
venient to use a perfect fluid analogy. The role of the 
pressure is played by the iagrangian p itself, while the 
energy density is given bya 



e = K(v)(2Xp, x (X) 

EE K^)i(X), 



P(X)) 



(4) 
(5) 



where ... t x denotes a partial derivative with respect to 
X. The ratio of pressure to energy density, which we 
call, for brevity, the fc-essence equation-of-state, 
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(6) 



does not depend on the function K(ip). For a "standard" 
kinetic term, p = X, in the case when there is no poten- 
tial, the equation-of-state is = 1. However, for a gen- 
eral choice of p it is easy to get any value of Wk ■ Notice 
that Wk < — 1 does not imply necessarily the instability of 
the fluid with respect to small wavelength perturbations. 
The effective "speed of sound," cs, which determines the 
propagation of perturbations in the fc-essence component 
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P,x 



P,x 



(7) 



and it can be positive for any Wk- For instance, the effec- 
tive speed of sound, defined to be the coefficient of the 
momentum-squared term in the perturbation equation 
for the scalar field, is always equal one for quintessence 
models with canonical kinetic energy, while the equation 
of state w can be rather arbitrary here. 

We want to study the evolution of a universe filled by 
fc-essence (labeled in the eqs. below by "k") and matter- 
radiation (labeled by "m" in cases where we refer generi- 
cally to the dominant matter-radiation component, be it 
dust-like or radiation, or by "d" or "r" if we refer specif- 
ically to the dust-like or radiation component, respec- 
tively). There is increasing evidence that the total encrgjz. 
density of the universe is equal to the critical value Jd'tEl 
and, hence, we will consider a flat universe only. In that 
case the equation for the scale factor a takes the form 



H 2 ee N 2 = e r 



(8) 



where a dot denotes derivative with respect to physi- 
cal time t and we introduced the number of e-foldings 
N = log a. This equation has to be supplemented with 
equations for e m and . These are the energy conserva- 
tion equations for each component j: 



(9) 



where Wj is the equation-of-state for the appropriate 
matter-radiation or fc-essence component. Considering 
a homogeneous field (p and substituting the expression 
for the energy density (^) into the appropriate Eq. (||) 
for fc-essence, one gets 



dX 
dN 



c 



3(l + w fe ) 



K, 



K H 



(10) 



where Wk is given by Eq. (g), a = sign(dtp / dN) and the 
Hubble constant is given by Eq. (||). 

We will consider functions p{X) that increase mono- 
tonically with X. They should satisfy further restric- 
tions, which follow from the requirements of positivity of 
the energy density, 



I = 2Xp :X - P > 



(11) 



and stability of the fc-essence background, c 2 s > 0, imply- 
ing 



e,x = 2Xp,xx +P,x > 0. 



(12) 



For designing models and visualizing constraints, it is 
helpful to re-express p as p = g (y) /y and consider it as 
a function of the new variable y = X^ 1 / 2 . The pressure 
of the fc-essence component is, therefore, 



p = K((f)g(y)/y; 



(13) 



the equation-of-state and the effective sound speed are, 
correspondingly, 
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(14) 



and the restrictions Eq. and (|T^) take the very simple 
form 



-g' > 0, e x = ^g" > 0, 



(15) 



where prime denotes derivative with respect to y. These 
conditions just mean that g should be a decreasing convex 
function of y — X~ x l 2 . A generic function which satisfies 
these restriction is shown in Fig. 1. Taking into account 
that H = 



./£tot = V £ m + £ k and e k = K (tp) e (y) = 
—Kg' (y) one can rewrite Eq. ( |To| ) in terms of the new 
variables as 



dy = 3 {w k (y) - f ) 
dN : 



r' (y) 



■(y) 



2K 3 / 2 



£k 
Stot 



where 

r (y) = (- 



V) 1/2 y (i + O 



3 (g-g'y) 
2V2 



(16) 



(17) 



is a function which, as we will see later, is critical for the 
attractor properties of fc— essence. 



3 




FIG. 1. A sample function g(y). Boldface letters denote 
the corresponding attractors; their positions are given on the 
t/-axis. The tangent to the curve at a radiation tracker, such 
as R, goes through 4j/_r/3, whereas the tangent to the curve 
at the de Sitter point S goes through the origin. 



III. CLASSIFICATION OF TRACKER AND 
ATTRACTOR SOLUTIONS 

The attractor solutions for fc-essence can be divided 
into two classes. In one class, fc-essence mimics the 
equation-of-state of the matter-radiation component in 
the universe. We refer to these as trackers because the 
cosmic evolution of fc-essence follows the track of an- 
other energy component. The second class of attractors 
consists of cases where fc-essence is drawn towards an 
equation-of-state which is different from matter or radia- 
tion. These attractors are important in the limits where 
fc-essence is either a negligibly small or an overwhelm- 
ing large fraction of the total energy density. The types 
of attractors available at any given moment in cosmic 
history depend on whether the universe is radiation- or 
matter-dominated. For all types of attractors, there is an 
associated basin of attraction, a set of initial conditions 
which evolve towards the attractor. 

In the presence of a matter background (dust or ra- 
diation) component with constant equation-of-state w m , 
Eq. ([Hi]) can have tracking solutions for which the fc- 
essence equation-of-state equals w m . To reveal when it 
can happen and to find these solutions explicitly we just 
need to note that if such solutions exist, they have to be 
generically of the form y (TV) 

— Dm — const., where y m 

satisfies the equation 



Wk (j/m) = — — ; 

yg 



(18) 



Substituting this ansatz into Eq. ( fig ) and noting that 
the ratio e k /et t should stay constant during the track- 
ing stage, we see that y (N) = y m can be a solution of 



Eq. (pf), only if K (tp) — const./ 1 ip 2 and, therefore, for 
simplicity, we consider from now on only scalar fields 
with Lagrangian 



g{y) 

(p 2 y ' 



(19) 



It is worth noting that this kind of dependence on a 
scalar field occurs in the string tree-level effective action 
when expressed in the Einstein frame.EirEj In this case, 
Eq. (pf) simplifies to 



dy _ 3 (w k (y) - 1) 
dN ~ 2 r> (y) 



Etot 



(20) 



where we restrict ourselves to the most interesting case 
of positive a on the branch of positive p. To close the 
system of equations for the two unknown variables y and 
£fc/e*ot, we use the equation 



d(e k /e t ot) 
dN 



£tot 



£tot / ^ 



wk(y)), (21) 



which immediately follows from Eq. 
lution of Eq. (pi), then y (N) = y, 
Eqs. <M) and (21), provided 



, = const, 



is a so- 
satisfics 



(Vr, 



\£tot 



< 1, 



(22) 



where the inequality is simply the physical constraint 
that e k < £tot (assuming positive energy densities 
and e m ). If r (y m ) > 1, a tracker solution y(N) = y m 
is physically forbidden. 



A. When are trackers attractors? 

To find out when trackers are stable solutions with a 
non-trivial basin of attraction, we study the behavior of 
small deviations from the tracker solution. Substituting 
y (N) = y m +Sy and e k /etot (N) = (ek/stot) m +S(e k /e t ot) 
into Eqs. ( po| ) and (^l]) and linearizing, we obtain 



dSy _ 3 (w k (y m ) - 1) 
dN ~ 2 rL 



5(e k /e tc 



(23) 



dS {ek/etot) 
dN 



-3r 2 m {l-r 2 m )w k {y m )6y, (24) 



where the index "m" denotes evaluation of the appropri- 
ate quantities at the tracker point y m and {e k /e t ot) m has 
been replaced by r 2 {y m ) according to Eq. d22|). Differ- 
entiating Eq. ( p3| ) with respect to N and using Eq. j24|), 
one obtains the following closed equation for 8y: 



d 2 Sy 

In 2 



+ I (1 - Wm) ^ + ?(1 - r 2 m ){l + w m ){c 2 s - w m )Sy = 0. 



(25) 
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Here c| is the squared "speed of sound" of k- essence 
at the tracker point and we took into account that 
Wk {Vm) = w m . Eq. ( |25| ) is a second order differential 
equation with constant coefficients and has two expo- 
nential solutions, ft is easy to see that for \w m \ < 1 both 
solutions decay if 



C S > Wr, 



(26) 



Therefore, since c 2 s 



(.9 — .<?'?/) /d"y 2 > an y tracker can 
be easily made an attractor by arranging a small second 
derivative of g at the tracker point. 

As important examples, let us consider the two most 
interesting cases, namely, trackers in the presence of radi- 
ation (labeled "r" in the equations below) and cold mat- 
ter (labeled "D" for "dust"). 



B. Radiation trackers 



For radiation trackers, w m = w r = 1/3 and Eq. (|18[), 
which defines the location of the radiation trackers (y m = 
y r ), reduces to 

Vrg'iVr) = -3.g(y r )- (27) 
The ratio of the energy densities is given by 



gfc \ 



r 2 (y r ) = -2g'{y r )y 2 r 



(28) 



and radiation trackers exist only if at the points y r sat- 
isfying Eq. (p7|), r 2 (y r ) < 1. These trackers are sta- 
ble attractors only if g" (y r ) < — 4g'(y r ) jy r . Radiation 
trackers are always located in the region where g > 
(positive pressure), corresponding to y < yn in Fig. 1. 
For a given g(y), there can be more than one radiation 
tracker. For each of them, the geometrical way of find- 
ing the value of y corresponding to the tracker is given 
in Fig. 1. These trackers can have different values of 
r2 {Vr) = (£fc/£iot) r - Numerically, a likely range for r 2 (y) 
is 10 _1 — 10 -2 . This is also the range we wish to have 
in order that cosmic acceleration begin at roughly the 
present epoch. We label the radiation tracker with the 
desired value of r 2 (yn) as R, and a second possible ra- 
diation tracker with a different value of r 2 {y r ) (the one 
closest to yr>) as r(?) in Fig. 1. If r 2 (y r ) is much smaller 
than 1(T 2 , the energy density falls so much at the on- 
set of matter-domination (before it freezes at a constant 
value) that it would not yet have overtaken the matter 
density today. If r 2 {y r ) is much greater than 10 _1 , then 
the contribution of fc-essence to the total energy density 
would change the expansion rate in the early universe 
and adversely affect the predictions of primordial nucle- 
osynthesis. The current constraints-pn r 2 (y r ) from rut. 
cleosynthesis vary from 4 per centta to 20 per cent,t£l 
depending on how the observations are weighted. 



C. Dust trackers 

The fc-essence field can also track the dust {wd = 0) in 
the (cold) matter-dominated universe. Since the pressure 
is proportional to g{y) and is zero for dust, it must be 
that 



g(yi 



0, 



(29) 



at the dust attractor point, y — yo- An additional condi- 
tion for the existence of the dust tracker is that r (yo) < 1 
(see discussion following Eq. (p2])). In this case the ratio 
of energy densities at the dust tracker is given by 



(-) 

\£tot/ D 



(y D ) 



■g'(yD)y 2 D - 



(30) 



If a dust tracker exists then it is always an attractor, 
since the stability condition Eq. (|2^) just means here 
that the "speed of sound" of fc— essence should be posi- 
tive. Note, that for the monotonically decreasing convex 
functions g under consideration only a maximum of one 
dust attractor can exist (see Fig.l) since g has only one 
zero. It is very important to point out that one can easily 
avoid a dust tracker by considering functions g such that 
r 2 (y D ) = -lg'{yD)y 2 D > 1 at y D , 



D. De Sitter Attractors 

We have noted that fc-essence can have attractor solu- 
tions which are not trackers in that they do not mimic 
matter or radiation. These attractor solutions play an 
important role in two extreme cases, namely, when the 
energy density of matter or radiation is either much big- 
ger or much smaller than the energy density of fc-essence. 
In this subsection, we study the case when the back- 
ground is dominated by matter-radiation and fc-essence 
is an insignificant component, Ek <C £ m - In this case, if 
g{y) satisfies some simple properties, fc-essence has an 
attractor solution in which it behaves like a cosmological 
constant (wk — > — 1). We refer to this solution as the de 
Sitter attractor (labeled "S"). 

Our purpose is to construct models in which fc-essence 
has a positive pressure, radiation tracker solution (R) 
during the radiation-dominated phase and approaches a 
state with negative pressure shortly after the onset of the 
matter-dominated phase. At the very least, it is neces- 
sary that g{y) be positive for some range of y and nega- 
tive for another range since the pressure is proportional 
to g{y). This simple condition is generically sufficient 
to produce a de Sitter attractor solution: Since g' must 
be negative (the positive energy condition, Eq. ([l5|)), 
it follows that g must have a unique zero, j/d, the only 
dust attractor possible. Furthermore, g(y) is positive for 
y < yo, a range which must include the radiation tracker, 
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y = ur. For y > yu, the pressure (oc g) and, corre- 
spondingly, uik = —g/yg 1 are negative. From this obser- 
vation, combined with the stability condition [g" > 0; 
see Eq. (15)), it follows that the derivative of r(y) (see 
definition (Il7j)) 



g"y 



7 (Wk - 1) 



(31) 



must be negative for y > yu- Since r (y) is positive at 
y = yu and has a negative derivative for y > yu, gener- 
ically (provided r' does not approach zero too rapidly) 
r(y) should vanish at some point y = y$ > yn and then 
become negative. As immediately follows from the defi- 
nition of r (see (p"7j)), the equation-of-state of fc-essence 
at y = ys (point S in Fig.l) corresponds to a cosmolog- 
ical term: Wk{ys) = — 1- Hence, we see that de Sitter 
attractors exist for a very wide class of g{y) and are a 
generic feature of fc-essence models. 

In the absence of matter, y (N) = ys = const, is not 
a solution of the equations of motion. However, when 
matter strongly dominates over fc-essence (sk/£tot €1), 
there exists a solution in the vicinity of this point. (For- 
mally, in the limit s^/stot ~ > 0, y (N) — > ys is an exact 
solution of Eqs. ( pp| ) and (fjl|).) Setting w m — Wk = — 1 
in Eq. ( p3[ ) it can be also verified that this is a stable 
attractor. For finite, but very small ratio Sk/stot *C 1, 
the approximate solution, corresponding to w rs — 1, is 
located in the vicinity of ys and has the form: 



— (AO oc exp (3 (1 + w m ) N) 

Stat 



and 



y (N) « y s + 



(3 + w m ) r' (y s ) \s to t 



— (N) 



1/2 



(32) 



(33) 



As shown below, if at any moment of time Ek jstot lies 
below the basin of attraction of the tracker solutions, k- 
essence will be driven first to the de Sitter attractor and 
stay in its vicinity as long as Sk/stot is sufficiently small. 
We will utilize this property at the transition from the 
radiation- to the matter-dominated phase. 



universe 



from Eqs. (17) and 



1 + w k (yje) 



1 The equation-of-state can be easily obtained 



2^2 



1 



3 V^M 

and the scale factor is 



= const, 



(35) 



(36) 



If —g' k yl/2 > 1 the solution describes power law infla- 
tion, which is an attractor of the system provided that 
r' (y k ) < 0. 

The existence of a fc-attractor depends mainly on the 
form of the function r(y). A /c-attractor corresponds to 
r {Vk) ~^ 1 (i.e., the limit where the energy density is 
totally dominated by fc-essence). In general, if r(y$) > 1 
for some yo and there exists an S-attractor (r(ys) = 0), 
then there must exist a fc-attractor somewhere between 
them, yo < yk < ys, simply because r (y) is a continuous 
function. 

In particular, we are interested in the case where there 
is no dust attractor because r(yu) > 1, and yet there 
is a de Sitter attractor with r(ys) = 0. In this case, 
not only must there exist a fc-attractor at some yo < 
Vk < ys 1 but we know that it has negative pressure (since 
9{vk) < 0), is stable (since Wk < 1, see Eq. (|3l"l)) and is 
the unique fc-attractor with negative pressure (since r' is 
monotonically decreasing in this y— interval). 

Note also that this negative-pressure fc-attractor only 
exists if there is no dust tracker solution, that is, r(yo) > 
1. If there is a dust tracker, (r(yo) < 1), then, since 
r'(y) < for y > yu, there is no point y — yx > yD 
where r(y) = 1 and, hence, there is no fc-attractor at 

Vd <y < ys- 

It is possible to have other fc-attractors with positive 
pressure at y < yjj (the closest one to yu is denoted by 
k(?) in Fig.l), but they will prove to be irrelevant in our 



scenario. 



IV. COSMIC EVOLUTION AND ATTRACTOR 
SOLUTIONS 



E. fc-Attractors 

Whereas the de Sitter attractors are important when 
fc-essence is an insignificant contribution to the total en- 
ergy density, the fc-attractors arise when fc-essence is the 
dominant energy component. In the absence of matter 
(£k/£tot = 1), the function y (N) = y k = const, where y k 
satisfies the equation 

r(y k ) = l, (34) 

is a solution of Eq. (^0|) , while Eq. ( ^T| ) is satisfied iden- 
tically. This solution describes a power-law expanding 



Once all possible attractors for fc-essence have been 
identified, it is easy to understand the evolution of the 
fc-field as a voyage from one attractor solution to another 
as different phases of cosmic evolution proceed. For both 
the radiation- and matter-dominated phases, there are 
several possible configurations of relevant attractor so- 
lutions. In this section, we systematically classify the 
attractor configurations for each phase and their conse- 
quences for cosmic evolution. 



G 



A. Radiation-Domination 

We assume that g(y) has been chosen so that there 
exists an attractor solution (R) at y = yn such that 
r 2 (yn) = e k /stot is in the range one to ten percent, 
roughly equipartition conditions. This energy ratio leads 
most naturally to a matter-dominated epoch that lasts 
a few billion years and cosmic acceleration beginning at 
about the present epoch. Depending on the form of r 2 (y) , 
which is determined by g(y) in the Lagrangian, there 
will be additional attractors during the radiation epoch. 
Whether y is drawn to the correct attractor y^ depends 
on initial conditions and the other attractors. Ideally, 
we want y = yn to have the largest basin of attraction 
so that most initial conditions join onto the desired cos- 
mic track. The combination of cosmologically relevant 
attractors during the radiation-dominated phase can be 
one of three types: 

A r ) R, S and no other attractors at ys > y > Vr- 
This occurs only if the function r (y) decreases for yu. < 
y < ys- Conversely, if r (y) increases somewhere in the 
range y > yn then it inevitably leads to the appearance 
of an extra k and/or r attractor at y > yR. Let us prove 
it. 

If the function r (y) increases within some interval, it 
means that the derivative r' (y) is positive there. On the 
other hand, as it follows from (|3l|), r'(y) is positive only if 
w k > 1. Since w k (y R ) = 1/3 , w k (y s ) = -1 and w k (y) > 
1 somewhere in the interval yn < y < yo, there must be 
another point y within this interval, where w k (y) = 1/3. 
If r (y) < 1, this point is a radiation tracker different from 
R with a different value of r 2 (y). If r 2 (y) > 1, then y is 
not a tracker at all; but, since r (ys) = — 1, there must 
exist a point in the interval ys > Vk > y where r (y k ) = 1, 
which corresponds to a fc-attractor. That is, either there 
is a an extra radiation tracker or there is an extra k- 
attractor. 

For models of type A r where r 2 (y) is monotonically 
decreasing, a dust tracker solution with r(yr>) < r(yn) is 
inevitable and k— essence will be attracted immediately 
to it after matter-radiation equality, a situation we are 
trying to avoid in order to explain the present-day cos- 
mic acceleration. The model p(X) = — 1 + X falls in the 
above category; with a field redefinition, the action can 
be recast into the model of a field with canonical kinetic 
energy rolling down an exponential potential)] an exam- 
ple which is well-known to track in both the radiation- 
and matter-dominated epochs. 

B r ) R, S, K plus possibly other attractors at y < y D . 
This situation takes place when there is no dust tracker 
solution (r (yo) > 1) the case considered in our first pa- 
per .0 

C r ) R, S (no K attractor) and at least one addi- 
tional attractors r(?) or k(?). This case occurs when- 
ever there is a dust tracker solution (r(yn) < 1)) with 



the property that r (yr>) > T (vr) or, in other words, 
(£k/£tot)D > {£k/£tot)R- Even though there exists a 
dust tracker solution, we will show it is nevertheless 
possible to have a finite period of cosmic acceleration 
at the present epoch before fc-essence reaches the dust 
tracker solution in the future. For this to occur, the 
function r (y) must increase somewhere in the interval 
Vr < V < Ud- This is precisely the case considered above 
(see discussion of case A r ), where we argued that there 
must be an extra r and/or fc-attractor in the interval 
Ud < V < Vr- Furthermore, the attractor closest to yo 
must have r(y r / k ) > r(yp) > r(yn); otherwise, we could 
find another attractor in the interval y r / k < y < yk, as 
can be shown by repeating the argument presented un- 
der A r for this interval. If r(y r / k ) > r(yu) > r(yn) 7 this 
second tracker has a larger fraction of fc-essence. 




FIG. 2. Phase diagram for case A r during the radia- 
tion-dominated epoch. Phase lines flow in the direction shown 
by the arrows, dashed horizontal lines determine the y coor- 
dinate of attractor solutions and boldface labels the corre- 
sponding attractor points. The dotted line shows the points 
where e k /e to t = r 2 (y). 

A phase diagram of the system of Eqs. (pO|)-(pl|) de- 
scribing the global evolution of the fc-ficld during radia- 
tion domination is shown in Figs. |2| || and |] for each of 
the cases A r , B r and C r respectively. Phase trajectories 
cannot cross the lines where £ k /stot is equal to zero or 
one, and, hence, their tangents are horizontal there. The 
position of the radiation tracker R is fixed by the intersec- 
tion of the y — yR line (dashed) and the r 2 (y) curve (dot- 
ted). If r 2 (y) is bigger than one at the intersection point, 
the tracker does not exist. Notice that the phase trajecto- 
ries go in the direction of increasing (decreasing) e k /etot 
for w k (y) < 1/3 (w k (y) > 1/3) and therefore, their tan- 
gents are horizontal at the points where w k (y) = 1/3. 
On the other hand, phase trajectories evolve in the di- 
rection of increasing (decreasing) y for s k /stot < f 2 (y) 
(ek/stot > f 2 {y)) and at the points where these phase 
lines cross the curve r 2 (y) their tangents are horizontal 
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(see Eq. (10)). The form of r(y) also gives a clue about 
the equation of state Wk(y)' in the region where r (y) is 
an increasing function of y we have Wk{y) > 1 and where 
it decreases Wk(y) < 1. Hence, as noted previously, r(y) 
is what mainly determines the structure of the phase di- 
agram. 

As clearly seen in the figures in all cases, if the fe-field 
is initially located near the R-trackcr, it converges to it. 
Therefore, the basin of attraction is non-zero in all three 
cases. The attraction region includes equipartition initial 
conditions, the most natural possibility. 

For A r , Fig. ||, the R-attractor has the largest basin 
of attraction, the complete phase plane. If one starts, for 
instance, at {e k /e to t)i = exp (-30) (efc/etot)fi, then the 
fe-field rapidly reaches the vicinity of the de Sitter point 
S and joins the attractor connecting this point to the 
R-tracker. 




FIG. 3. Phase diagram of a model of the type B r during 
the radiation-dominated phase. In the relevant region of the 
diagram all trajectories can be traced back to a common ori- 
gin. Some of the phase trajectories converge to the radiation 
tracker R, while others, after approaching the de Sitter point 
S finally reach the K-attractor. The saddle point x "sepa- 
rates" both types of trajectories. 

The cases B r and C r have limited basins of attraction, 
and so are not as favorable from the point of view of ini- 
tial conditions. If the energy density of the fe-field is much 
smaller than the value at the R-tracker, the fe-field trav- 
els first to the vicinity of the S-attractor, where it meets 
the phase trajectory that connects it to the K-attractor 
(case B r ) or the r-attractor (case C r ). In either situation, 
the field never reaches the R-tracker. Although the lat- 
ter two cases have smaller basins of attraction than case 
A r , only cases B r and C r can produce cosmic accelera- 
tion today. One can simply assume that the initial value 
of the fe-field lies somewhere in the basin of attraction, 
a reasonable possibility. An alternative is to introduce 
additional ^—dependence in the Lagrangian, as for in- 
stance, L = g(y, ip) /yip 2 , where g (y, <p) -> g x (y) at high 
energies (ip is smaller than some ipo) and g (y, ip) — > g-i (y) 



at relatively low energies (ip is bigger than ipo), such that 
<7i (y) has an A r -set of attractors and g2 (y) has a B r /C r - 
set of attractors. Note, that the exact value of <po is 
not important at all, we only have to be sure that the 
transition from one regime to the other happens before 
equipartition. Although modifying the Lagrangian may 
seem more complicated, it has the advantage that it re- 
moves nearly altogether dependence on initial conditions. 



r 2 (y) 




FIG. 4. Phase diagram of a model of the type C r during 
radiation domination, with same notation as in Fig. 3. 



B. Matter Domination 

We have shown that it is possible to choose a wide 
range of models and initial conditions for which the fe- 
field converges to the R-tracker during the radiation- 
dominated epoch. The goal is to produce a scenario in 
which fc-essence overtakes the matter density and induces 
cosmic acceleration today. Yet, the contribution of fc- 
essence to the total energy density must not spoil big 
bang nucleosynthesis or dominate over the matter den- 
sity at the end of the radiation-dominated epoch (see Sec. 
III.B). To satisfy these conditions, it typically suffices if 
the R-tracker satisfies 



{ek/et 



ot)R 



r 2 (vr) 



a 



10' 



10" 



(37) 



In this subsection, we study the evolution as the uni- 
verse enters the matter-dominated epoch and the fc-field 
is forced to leave the radiation tracker. In a dust dom- 
inated epoch the relevant attractors can appear in the 
following two possible sets: A^) S, K and B^) S, D. 

In both cases successful fc-essence models are possible. 
In the case Ad, which was discussed in our earlier paper ,B 
there is no dust tracker solution, (r (yn) > !)• There- 
fore, as seen in the phase diagram of Fig. ^, when the 
radiation-dominated epoch is over, fc-essence approaches 
first the S-attractor; afterwards, when its energy density 
has increased significantly, it moves to the K- attrac- 
tor (a state with negative pressure but Wk > — !)■ If 
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Wk (uk) < —1/3, the expansion rate accelerates for the 
indefinite future; if —1/3 < Wk [uk) < 0, the expansion 
rate decelerates. Either way, the matter-radiation den- 
sity is increasingly negligible compared to fc-essence in 
the far future. 

In the second case (B^), there is a dust tracker so- 
lution. If (sk/£tot)D <C 1, fc-essence contributes only a 
small fraction of the total energy density at this attractor, 
and it approaches this attractor almost immediately af- 
ter matter-radiation equality. This is not desirable since 
then fc-essence cannot dominate today or cause cosmic 
acceleration. However, if (sk/etot)D = r 2 (yo) — > 1 or 
(efc/£d)r> ^> 1, there can be a period of cosmic accelera- 
tion before the fc-field reaches the dust attractor since it 
can first approach the S-attractor and remain there for 
a finite time, see Figure 6. Ultimately though, the ac- 
celeration is temporary; the fc-field proceeds to the dust 
tracker, the expansion of the universe begins to decel- 
erate, and the ordinary and (cold) dark matter density 
approaches a fixed, finite fraction of the total energy We 
refer to the scenario as a "late dust tracker" because the 
dust attractor is reached long after matter-domination 
has begun. 

Taking into account that r{yo) is near unity or greater 
for both case and B^, we obtain from Eqs. (17) and 

&■ 



9'rVr 

q'd v 2 d 



< - a~5-(l(r 3 -l(r 2 ). 
16 



We can also infer from Fig. 1 that g' D ■ (y^ 
9{vr) — ~Vr9r/^ and, therefore, for 



VR ^ 3 
— < — a 

yo 16 



2 • (1(T 3 - 1CT 2 ) 



and 



lR < " ~ 6 . (io-4 _ io-3). 
9r 16 



(38) 



Vd) < 



(39) 



(40) 



Since e k = -g'/<p 2 and \g' (y s )\ < \g' (vd)\, we con- 
clude that after radiation domination, when the fc-field 
reaches the vicinity of the S-attractor, the ratio of en- 
ergy densities in fc-essence and dust can not exceed 
Sk/sd < a 2 /l6 ~ 6 • (10- 6 - 10- 4 ). This is the nadir 
of fc-essence; once fc-essence approaches the S-attractor, 
its contribution to the cosmic density increases again un- 
til it becomes comparable to the matter density. In case 
Ad, the fc-field will evolve further to the K-attractor and 
the fc-essence energy will increasingly dominate over the 
matter density. In case B^, the fc-field approaches the D- 
tracker where the ratio of fc-essence to the matter density 
approaches some fixed positive value. 

The statements above are generic and do not depen- 
dent significantly on the concrete model as long as it 
satisfies the simple criteria formulated above. Let us 



stress that the only "small" parameter used is the ratio 
(e/c/£tot)_R, which has to be of the order of lO -2 — 10 _1 , a 
very natural range for these models and one that satisfies 
constraints of big bang nucleosynthesis (see Sec. III.B). 
For this range, the present moment is approximately the 
earliest possible time when cosmic acceleration could oc- 
cur. 




FIG. 5. Phase diagram of a model of type Ad during the 
matter-dominated epoch. All trajectories have a common ori- 
gin and all of them finally reach the K-tracker. Trajectories 
which "skim" the line Ek/etot ~ reach this attractor after 
going through a nearly de Sitter stage (the S-attractor). 




FIG. 6. Phase diagram of a model of type Bd during the 
matter-dominated epoch. All trajectories have a common ori- 
gin and all of them finally reach the D-tracker. Trajectories 
which "skim" the line Ek/etot ~ reach this attractor after 
going through a nearly de Sitter stage. 

Finally note that, during the transition from the radia- 
tion tracker R to the de Sitter attractor S, the equation- 
of-state of fc-essence has to take values bigger than one, 
and hence the dominant energy condition Sk > \pk I is vio- 
lated during a certain finite time interval. This violation 
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implies thaL-fc-essence energy can travel with superlumi- 
nal speeds.li-3 Thus, perfectly Lorentz-invariant theories 
containing non-standard kinetic terms seem to allow the 
presence of superluminal speeds, as already pointed out 
in.BB 

V. CONSTRUCTING MODELS 

In previous sections, we have presented a general the- 
oretical treatment of the attractor behavior of fc-essence 
fields in a cosmological background. We have empha- 
sized the properties needed to formulate models which 
will lead naturally to cosmic acceleration at the present 
epoch. In this section, we discuss how to apply the gen- 
eral principles to construct illustrative toy models. 

Let us summarize the conditions we have derived for 
building viable Lagrangians. First, we must satisfy 
the gen eral positive energy and stability conditions in 
Eq. (|l5|). If g takes positive and negative values, they 
already suffice to guarantee generically the existence of 
a radiation point y R where w(yn) = 1/3, a unique dust 
point yn where w(yo) = 0, and a unique de Sitter point 
ys where w(ys) = — 1. The radiation point is an attrac- 
tor if g"(yn) is sufficiently small, 

g » ( y R) < _ 4 ^M, (41) 

vr 

and the remaining prerequisites needed to ensure a suc- 
cessful scenario are then reduced to simple restrictions 
on the derivative of g at two separate values of y: 

i) At y R , r\ = -2g'(y R )y 2 R ~ 1CT 2 - 10" 1 . 

ii) At y D either r 2 D = -9 y 2 D g'{ijD)/& > 1 or 1 - r 2 D = 
i + 9yl 5'M/8<1. 

The first condition in ii) corresponds to cases where there 
is no dust attractor, and the second condition to cases 
where there is a dust attractor with a small matter to 
fc-essence energy density ratio. 

A straightforward way of constructing a function with 
given derivatives at two points is to glue two linear func- 
tions with the required slopes, as shown in Figure |7[ Ob- 
serve that if g{y) is linear around the radiation point the 
attractor requirement ([ll]) is automatically fulfilled. In 
order to have a finite c|, it suffices to introduce small 
quadratic corrections to the glued linear functions. We 
implement this procedure to build a toy model expressed 
in terms of artificial parameters (from the point-of-view 
of fundamental physics) that can be simply related to 
Fig. |?] and our earlier discussion of attractor solutions. 
One should appreciate that, for this pedagogical purpose, 
we have "overparameterized" the problem - the outcome 
is rather insensitive to most parameters as long as they 
obey certain simple general conditions. Simpler forms 
with fewer parameters are certainly possible. 



g(y) 



ii 




FIG. 7. A simple toy model for g(y) consisting of two 
linear pieces meeting at the "crossing point" y c . Here yn and 
yD are the radiation and the dust attractor values, and the 
derivatives of g at these points are g' R and g' D , respectively. 

Let g g iue{y) be any smooth function constructed by 
gluing the two linear pieces of Fig. [?]. The function g g i ue 
depends on y and has yR,g' R ,yD and g' D as parameters 
where y R and yo are the radiation and the dust attractor 
values and the derivatives of g at these points are g' R and 
g' D respectively. Our toy model corresponds to 

g(y)=g glue (y)(l-^y—]. (42) 

The factor g g i ue describes the function in Figure ^ and the 
factor in parenthesis provides the quadratic corrections 
needed to have a positive speed of sound. It so happens 
that the latter factor also shifts the de Sitter point from 
y = oo, as it would be for purely linear functions, to finite 
y, although this is not crucial for our purpose. For s>l 
the de Sitter point is located at ys ~ s-yo and g « g g i ue . 

Once a general form for g is known, such as the ex- 
ample above, one can study how the model parameters 
affect the resulting cosmology. Our conclusion is that 
the predictions of the toy model are relatively insensi- 
tive to the gluing function or to the particular values of 
ZARj g'n-i g'n an d ys as long as they satisfy certain sim- 
ple relations. For instance, what sets the values of 
and Wk today? Do these depend on the precise form of 
the interpolating function? We have solved numerically 
the equations-of-motion for a wide range of gluing func- 
tions g g iue in Eq. (0). For a typical parameter choice, 
the final value of £lk does not depend on the particular 
gluing function as long as g g i ue conforms closely enough 
to Fig. 

The value of today does depend on the evolution of 
Sk/s m . At early times the field is locked at the radiation 
tracker, and its fractional energy density ratio is given 
by —2g' R y 2 R . After radiation-matter equality the field can 
not follow the radiation tracker anymore and its energy 
density drops by several orders of magnitude until et / e m 
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reaches a minimum value at the time Wk falls below zero. 
We shall label this minimum value with the subscript 
"min" . The energy density at this minimum is roughly 
given by 

(?) * (43) 

The position of the minimum in time only depends on 
the distance between the radiation and crossing point 
Vc — Ur- As y c — y R increases from zero, the minimum 
is shifted from matter-radiation equality to later times. 
After reaching the minimum, the field moves on to the 
de Sitter attractor and Ek/e m grows as (z + 1)~ 3 , where 
z is the red shift. In order to have fc-essence dominate 
today, it must be that £k / e m during the radiation epoch 
lies roughly between 10 _1 and 10~ 2 . Then, (sk/£ m )min 
lies in the range 10~ 4 — 10~ 6 and, provided y c is chosen 
appropriately this has fc-essence dominating at about the 
present epoch. One can see these conditions impose con- 
straints on certain combinations of our parameters, al- 
though in a fairly natural range not very far from unity. 

As discussed in Section IV. B, there are two possible 
future fates for the universe depending upon whether 
there is a "late dust tracker" solution or not. By re- 
quiring r 2 D > 1 we avoid a dust tracker and, therefore, 
insure that the fc-field approaches the fc-attractor when 
fc-essence starts to dominate. The equation of state of 
fc-essence at the fc-attractor depends on the parameter s. 
By increasing s the equation of state Wk at the fc-attractor 
approaches —1, and in the limit s — > oo, Wk(jjK) ~~ * — 1- 
If Wk < —1/3, the expansion rate of the universe ac- 
celerates forever. Using the maximal value of the w at 
the present epoch as allowed by supernovae observations, 
say, s can be simply adjusted to insure that w at the fc- 
attractor is less than or comparable to this value. In this 
case, the equation-of-state of fc-essence today will be less 
than or equal to Wkiyx), which is set by s, as described 
above. 

If Tjj < 1 , it is possible to have successful models if r 2 D 
is sufficiently close to 1. In such a model the equation- 
of-state of fc-essence will finally reach Wk = in the far 
future; so, ultimately, cosmic acceleration ceases and the 
expansion begins to decelerate again. Nevertheless, it is 
still possible to have a finite period in which the equation- 
of-state is negative and which includes the present epoch. 
It is worth noting that models without a dust attractor 
are more generic and natural, since they do not require a 
special tuning of r (j/d) to a value close but smaller than 
unity at the dust point. Below we illustrate examples of 
both types. 



A. Model without dust attractor 
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FIG. 8. The ratio of fc-essence to matter energy density, 
Sk/£m, vs. 1 + z for a model with a fc-attractor. 
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FIG. 9. The equation-of-state Wk vs. 1 + z for a model 
with a fc-attractor. 

Models that belong to the general class A^ illustrated 
in Fig. 5 do not have dust attractor solutions because 
K2Ad) > 1- Choosing the following values of the param- 
eters, y R = 0.1, g' R = -5, y D = 17, g' D = -5 • 10~ 3 and 
s 2 ■ yo — 135, we have r(yo) ~ 1.2. Therefore, there 
has to be a K-inflationary attractor, which is located for 
our parameter choice at yx ~ 28. At the K-attractor, fc- 
essence has the equation-of-state Wk {vk) — —0.43. The 
ratio of the energy densities at the R-tracker in this 
model is (sk/ £ tot)R = 0.1. The results of the numeri- 
cal calculations are presented in Figs. || and ||. We see 
that during the radiation stage fc-essence quickly reaches 
the radiation tracker, in particular, the oscillations of the 
equation-of-state Wk in Fig. ^ around Wk — 1/3 decay 
exponentially rapidly. The fc-field has the same equa- 
tion of state as radiation until the moment when dust 
starts to dominate. Around this time the energy density 
of fc-essence suddenly drops by three orders of magni- 
tude and the equation-of-state, after a very short period 
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of increase, drops down to Wk — — 1, the value of the 
equation-of-state along the S-attractor. After that, when 
the energy density of fc-essence becomes significant, Wk 
starts to increase towards the K-attractor value, -0.43. 
Since fife is not yet unity, the current value is somewhere 
between the K-attractor value and -1; in this example, 
the value today (z = 0) is Wk — —0.69. The energy 
density of fc-essence today is fife « 0.65, and because we 
assumed a flat universe, fi m = 0.35. For completeness 
let us mention that we have defined "today" (z = 0) to 
be the moment when the matter-radiation energy density 
ratio is given by (e r /e m ) t oday = 4.307 ■ 10~ 5 /(fi m h, 2 ). 

B. Model with a late dust attractor 
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FIG. 10. The ratio of fc-essence to matter energy density, 
£fc/£m, vs. 1 + z for a model with a late dust tracker solution. 
In this type of model, Wk — * in the far future and the ratio 
of fc-essence to matter energy density approaches a constant. 
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FIG. 11. The equation-of-state Wj, vs. 1 + z for a model 
with a late dust tracker solution. 

Taking y R = 11 • 10" 3 , g' R = -34, y D = 11, g' D = 
—8 • 10 -3 and s 2 - yn = 56, we can construct a model with 
a "late dust tracker" , corresponding to the phase diagram 
in Fig. 6. The parameters have been deliberately chosen 



to differ significantly from the ones in the model without 
dust attractor in order to illustrate that fine tuning is 
not necessary. 

The late dust attractor is reached after fc-essence 
passes near the de Sitter attractor following matter- 
radiation equality. At the late dust tracker {Ek/£tot)D — 
r 2 (y D ) ~ 0.88 and, correspondingly, {e k /e d )D ^ 7. 
Hence, the fractional contribution of the matter density is 
small but remains finite in the indefinite future. The ra- 
tio of energies at the R-tracker is {ek/etotjR — 8-3- 10~ 3 . 
The results of the numerical calculations are presented in 
Figs. and[Tl|. The evolution of the fc-field here is very 
similar to the one we described in the previous case; the 
differences between both models occur at small red-shifts. 
The fraction of the critical energy density of fc-essence 
today is in this model also fife = 0.65 and the equation- 
of-state Wk takes the value —0.4. The future evolution 
of the model with a late dust attractor is completely dif- 
ferent from what we found in the previous one. Here the 
ratio of the energy densities of fc-essence and dust will 
continue growing in the future only until it becomes ap- 
proximately 7. After that it will start to oscillate around 
this value with exponentially decaying amplitude while 
the pressure approaches the dust point, where Wk = 0. 

C. Simpler and More Practical Examples 

The toy models presented thus far are all built on the 
ansatz shown in Fig. [?], which entails numerous parame- 
ters. We have pointed out that the large number of pa- 
rameters is not a necessary feature. We have introduced 
this form for pedagogical purposes, since it enables one 
to study directly the relation between the attractor solu- 
tions and cosmic evolution. Indeed, our analysis showed 
that the cosmological solution is relatively insensitive to 
most of the parameters provided they obey a few broad 
conditions. 

To emphasize the point, consider a model of the form 

p(X) = -b+ 2^1 + X -h(aX), (44) 

where h(aX) is some smooth function that can be ex- 
panded in a power series in X. This particular form 
is reminiscent of a Born-Infeld action in which h(aX) 
could represent higher order corrections in X. (This 
choice of a square-root form is not essential - simply 
an example.) As a specific case, for b = —2.05 and 
X-h(aX) = X-{aX) 2 + (aXf-{aX) A +{aXf-{aX/2f 
the Lagrangian defined by (0) satisfies all constraints 
and produces fi m = 0.3 and Wk = —0.8 today if one 
chooses a = 10~ 4 . This particular example has a cosmic 
evolution similar to the one described in Section V.A (no 
dust attractor). We see that in this case, as with a wide 
range of other functional forms, the condition b > 2 and 
the choice of the single parameter a suffices to satisfy all 
of the conditions of the multi-parameter toy models. 
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VI. DISCUSSION 

Introducing a dark energy component with negative 
pressure has resolved many observational problems with 
the standard cold dark matter model including the recent 
evidence from supernovae searches that the universe is 
undergoing cosmic acceleration. At the same time, the 
dark energy component presents a profound challenge to 
cosmology and fundamental physics. What is its compo- 
sition and why has it become an important contribution 
to the energy density of the universe only recently? 

The example of fc-essence shows that it is possible to 
find a predictive, dynamical explanation that does not 
rely on coincidence or the anthropic principle. Unlike a 
cosmological constant or quintessence models of the past, 
the energy density today is not fixed by finely-tuning 
the vacuum density or other model parameters. Rather, 
the energy density today is forced to be comparable to 
the matter density today because of the dynamical inter- 
action between the fc-essence field and the cosmological 
background. 

Technically, the fc-essence approach, at least in the ex- 
amples we have constructed, relies on attractor properties 
that naturally arise if the action contains terms that de- 
pend non-linearly on the gradients of the fc-essence field. 
Non-linear terms of this type appear in most models uni- 
fying gravity with other particle forces, including super- 
gravity and superstring models. In the past, these con- 
tributions have been ignored for reasons of "simplicity." . 
The example of fc-essence demonstrates that the effects of 
non-linear dynamics can be dramatic. In a cosmological 
setting, we have shown how they can cause the fc-essence 
field to transform from a tracking background field during 
a radiation-dominated epoch into a an effective cosmo- 
logical constant at the onset of matter-domination. This 
effect explains naturally why cosmic acceleration could 
begin only at low temperatures, at roughly the present 
epoch. 

The non-linear dynamics is totally missed if the kinetic 
energy terms are truncated at the lowest order contri- 
butions. Hence, the kinds of attractor effects discussed 
in this paper have gone unnoticed in most treatments of 
quantum field theory. This was one of the reasons for pro- 
viding a detailed, pedagogical treatment for at least one 
class of models. Clearly, this is the tip of a broad arena 
of study. As another possible application, it is interesting 
to note that a fundamental problem of superstring mod- 
els is to control the behavior of the many moduli fields 
in the theory, which are coupled to one another through 
non-linear kinetic energy terms. At the linear level, the 
moduli appear to be free fields with a flat potential, and 
so there is no guidance as to why, amongst all the possi- 
ble limits of M-theory, the low energy limit looks like the 
Standard Model. Perhaps non-linear attractor behavior 
constrains the evolution of moduli fields. 



In this paper, we have focused on how non-linear dy- 
namics addresses a fundamental theoretical issue, the 
cosmic coincidence problem. An important question to 
consider is whether there are observational tests to dis- 
tinguish fc-essence from alternative explanations. One 
notable feature of fc-essence models compared to the 
more general tracker quintessence modelsM is that the 
equation-of-state, Wk is increasing at the present epoch. 
For quintessence scalar fields rolling down tracker poten- 
tials, the quintessence tracks the matter density (w = 0) 
during most of the matter-dominated epoch, and only 
recently has begun to decrease towards w = — 1. Hence, 
measurements of dw/dz for the dark energy would distin- 
guish these two possibilities from one another and from a 
cosmological constant. However, this test would not dis- 
tinguish fc-essence from more general quintessence mod- 
els that can also be tuned so that Wk is increasing today 
as well. A second feature of fc-essence is the non-linear 
kinetic energy contribution. A consequence is that the ef- 
fective sound speed c 2 s is generically different from unity, 
whereas c s = 1 for a scalar field rolling down a potential. 
Depending on the model, the distinctive sound speed can 
have subtle or significant effects on the cosmic microwave 
background anisotropy. We will address these observa- 
tional considerations in a forthcoming paper J13 

As regards the future of the universe, our work here 
offers a new, perhaps pleasing possibility. In previous 
models with cosmological constant or quintessence, the 
acceleration of the universe continues forever and ordi- 
nary matter that composes stars, planets and life as we 
know it becomes a rapidly shrinking fraction of the en- 
ergy density of the universe. In the "late dust tracker" 
scenario which we have introduced here, the acceleration 
is temporary and the matter density approaches a fixed, 
finite fraction of the total. 
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